Dynamic image potential at an Al(lll) surface 
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We evaluate the electronic self-energy at an Al(lll) surface using the GW space-time method. 

This self-energy automatically includes the image potential Vim not present in any local-density 
approximation for exchange and correlation. We solve the energy-dependent quasiparticle equations 
to obtain surface state wavef unctions, and calculate the effective local potential experienced by 
electrons in the near-surface region. We find that Vim for unoccupied states is due to correlation (not 
exchange). The image-plane position for interacting electrons is considerably closer to the surface 
than for the purely electrostatic effects felt by test charges, and, like its classical counterpart, is 
drawn inwards by the effects of atomic structure. 



Electrons outside a metal surface experience a sur- 
face barrier which has the asymptotic form of an im- 
age potential Vim — —1/4(2 — Zq) where zq is the 
effective edge of the metal. For interacting electrons 
this is an exchange-correlation (XC) effect, the quantum- 
mechanical analogue of the charge density redistribution 
which gives the classical image force. The form of the 
surface barrier is important for interpretation of low- 
energy electron diffraction |Q and scanning tunneling mi- 
croscopy Q experiments. In addition, surface and image 
states bound by Vim can be directly observed by mod- 
ern inverse photoemission and two-photon photoemission 
experiments Q]. 

The physics of the electron-surface interaction, and 
its transition from quantum-mechanical behaviour to the 
classical limit, has been the subject of a wealth of theo- 
retical studies. It is well known that the effective poten- 
tial of density-fimctional theory (DFT) within the local 
density approximation (LDA) (or gradient-corrected ver- 
sions) fails to reproduce the image tail shape which 
is the result of long-range many-body effects. The clas- 
sical response of metal surfaces to an electric field or 
distinguishable test charges has been investigated us- 
ing DFT but the resulting image-plane position 
Zc may differ from that experienced by electrons. The 
quantum-mechanical XC potential has been investigated 
by Eguiluz et al. for a jellium surface while for real 
materials it has been shown that atomic structure beyond 
the jellium model plays an important role in determin- 
ing both Zc and the binding energy of image states | [Tc| ] 
(whose existence requires the presence of a surface band- 
gap not exhibited by jellium). 

In this Letter we present the results of a calculation 
of the non-local electronic self-energy 'S{r,r',E) for an 
Al(lll) surface, including its full variation with energy, 
evaluated within the GW approximation [ pd| , which al- 
lows treatment of long-range correlation effects from first 
principles. This work is the first application of the GW 
space-time method, outlined in a recent Letter fl^ , to a 
system with a large supercell. In general surface states 



are poorly described in the LDA (and image states and 
resonances are entirely absent) and we therefore solve 
the quasiparticle (QP) equations self-consistently for one- 
electron-like excitation energies, without recourse to first- 
order perturbation theory. The effect of E is interpreted 
in terms of the effective local potential Vioc felt by QP 
states in a given energy range. Vioc automatically con- 
tains the image potential and the form of the crossover 
from image to bulk behaviour which is crucial for states 
localised near the surface. We find zq (for quantum- 
mechanical electrons) at Al(lll) to be significantly closer 
to the surface than Zc (for external test charges). For 
Al(lll) Zq is shifted inwards relative to jellium by the 
atomic nature of the surface (as is seen for Zc). There has 
been ongoing controversy in the literature over the sepa- 
rate contributions of exchange and correlation 
and we show that the image potential felt by excited-state 
electrons is not present in an exchange-only calculation. 

In many-body theory, exchange and correlation are de- 
scribed by the self-energy 'S{r,r',E), which is non-local 
and therefore in principle state-dependent, varies with 
energy, and is complex, containing information about the 
lifetimes of QP excitations. The non-local effects which 
give the image potential are present implicitly in the 
exact exchange-correlation potential V^c of Kohn-Sham 
DFT (as opposed to V^J^"^), but the surface barrier felt 
by excited states and the lifetimes of surface states are 
two features which cannot formally be addressed within 
DFT. The GW approximation for the self-energy, written 
in real space and time, is 

I](r, r', t) = iG{r, r', t)W{r, r', t) (1) 

where G is the one-particle Green's function, and W 
is the dynamically screened Coulomb interaction. This 
first-order diagram for S has been shown to work well 
in a wide variety of studies of real materials, successfully 
predicting, for example, the band gaps of semiconductors 
and insulators and the valence bandwidths of simple 
metals |0 . We follow the usual procedure of performing 
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an LDA-DFT calculation to calculate a non-interacting 
approximation for G, and obtain W within the random- 
phase approximation (RPA). We note that W calculated 
at this level already contains image-like interactions, as 
RPA screening corresponds to time-dependent Hartree 
theory. 

In the space-time method S is constructed in real space 
and (imaginary) time, which is advantageous because the 
GW self-energy is then a product rather than a con- 
volution in reciprocal space and energy. In imaginary 
time, the structure of the many-body response functions 
is much smoother and thus well suited to numerical work. 
Fast Fourier transforms (FFTs) are used extensively to 
move between real and reciprocal space, and between 
imaginary time and energy. Our calculation begins with 
the formation of the Green's function G(r, r',iT) in real 
space and imaginary time. We then proceed via the non- 
interacting density response function and dielectric ma- 
trix to calculate W without the use of any plasmon-pole 
approximation for frequency dependence. E(r, r',«T) is 
then formed and matrix elements in the LDA eigenfunc- 
tion basis are computed, Fourier transformed to imagi- 
nary energy, and fitted to the multiple-pole form 

(V^„k|SMI^„'k) = aL'k + E (2) 

This form facilitates analytic continuation to the real en- 
ergy axis, and represents a highly controlled approxima- 
tion since the accuracy of the fitting of the calculated 
self-energy can be directly monitored. 

The system we study is an Al(lll) surface. We employ 
a slab geometry with 5 layers of aluminium (sufficient to 
give, for example, a well-converged surface energy and 
Zc) and 8 layers of vacuum. The LDA calculation for 
the slab was undertaken with a 4x4x1 k-point mesh and 
a high energy cutoff. In order to converge our surface 
barrier, the required parameters for the GW calculation 
were found to be a plane- wave energy cutoff of 9 Ry (cor- 
responding to a 5 X 5 X 60 real space grid in the unit cell) , 
and 243 imaginary time points with At = tt/IO a.u. 300 
bands (those up to an energy of 74 eV above Ep in the 
LDA calculation) were included in G. 

QP energies for bulk materials are usually evaluated 
within first-order perturbation theory, employing the as- 
sumption that the wavcfunctions given by the solution of 
the QP equations, 

(^-^V2 + 14xt(r) + V„{r)^ *„k(r) 

+ J dr'l](r,r',^„k)*nk(r') = ^„k*„k(r), (3) 

are sufficiently similar to their DFT Kohn-Sham (KS) 
counterparts, where S is replaced by Vxc- The ability 
to solve the full QP equation for real materials is poten- 
tially important for highly inhomogeneous systems such 



as heterostructures or defects where the LDA may give 
qualitatively incorrect states. 

In the case of the metal surface, the QP eigenfunc- 
tions will include surface and image states, bound by 
Vim, which will differ significantly from the LDA states. 
We therefore obtained QP eigenfunctions by diagonalis- 
ing the QP Hamiltonian in the LDA eigenfunction ba- 
sis at a trial energy, and then iterating the energy of 
the given QP state to self-consistency. The full energy- 
dependence of the self-energy matrix in the basis of LDA 
states is therefore required, making use of the space-time 
approach especially important, as a functional form for 
the energy dependence of S is found. Even for the very 
inhomogeneous surface-slab system, two poles give a sta- 
ble fit of extremely high quality (rms error 0.2%) for the 
diagonal matrix elements. Most off-diagonal matrix ele- 
ments are zero by symmetry, with the (sparse) remainder 
similarly well described by the two-pole form. In Fig. 
|l| we show an unoccupied surface state below the vac- 
uum level at F, obtained directly from solution of the 
QP equation. The weight in the near-surface region is 
significantly enhanced relative to the LDA state by the 
improved description of XC effects. 

Although S is non-local it can also be viewed as a state- 
dependent local potential. In particular, in the asymp- 
totic limit far from the surface we expect the effect of 
S to be that of an image potential, independent of QP 
energy ^7^. Comparing the Hamiltonian for a local and 
non-local potential, it is clear that the state-dependent 
effective local XC potential Vioc for a QP state is defined 
by 

Vioc{r)^Qp{r)^ J dv'j:ir,r\EQp)^Qp{v'), (4) 

as used by Deisz et al. jisj in their study of a jellium sur- 
face. In the space-time method S is obtained on the real 
energy axis in the form of matrix elements in the LDA 
eigenfunction basis. We therefore use the completeness 
relation to write 

Vioc{v) - E<^I^")(^"I^(^Q^)I*QP)/<'"I*QP)' (5) 

n 

where the sum is over all KS states tpn at a given k-point. 
As the potential is not defined by Eq. || at the nodes of 
the QP state, we take a weighted average of the resulting 
V;oc(r), according to \9qp{v)\^, over a few states in a 
small energy range. 

The resulting surface barrier, split into contributions 
from exchange and correlation, is shown in Fig. |^. This 
potential was calculated using four states at the F point 
within 1 eV of the vacuum energy, and was well con- 
verged with the inclusion of 200 bands in the sum over 
off-diagonal matrix elements. 

As expected, the bulk value of Vjoc is similar to V^^^"*, 
as the QP energy shift for states near Ep in bulk 
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Al is rather small. Moving out through the surface, 
where V^/'^ falls exponentially (as the density does), 
Vioc crosses smoothly to the asymptotic image potential 
(shown with the best-fit zq). This image potential has 
been modified to take account of the repeated slab ge- 
ometry, which gives rise to two infinite series of image 
charges, but the resulting form becomes very similar to 
an isolated image potential within 10 a.u. of the surface. 
The form of the crossover is often treated in an ad hoc 
manner which is somewhat arbitrary for Al(lll) as the 
classical Vim and V^J^^ do not meet. However it is in- 
teresting to note that the dynamic image potential limit 
which we have calculated in this work comes much closer 
to meeting V^J^^. 

Surface corrugation of the XC potential, which is 
well described by the LDA near the surface, quickly 
decays outside the surface. Our results are consistent 
with a surface-position-independent image plane, as ar- 
gued and demonstrated for test charges by Finnis et 
al. [||. The potential also proves to be almost entirely 
state- independent; in the case of the image tail this en- 
ergy independence arises through cancelling effects from 
the spatial character of different states and the energy- 
dependence of S. 

The 'classical' image-plane position Zc has been cal- 
culated for jellium, and more recently for real materi- 
als, using a variety of techniques, generally based on the 
response of the ground-state charge density to external 
fields PjTHigj] and test charges ||]. The image-plane po- 
sition has also been found through the self-consistent re- 
sponse of a modified effective potential |^^. These re- 
sults, together with the value for dynamic electrons 
at a jellium surface obtained by Eguiluz et al. j9j, are 
compared in Table | with zq for electrons at an Al(lll) 
surface as calculated in this work, where the geometric 
edge is half a layer spacing outside the outermost crystal 
plane. 

Our value for zq is closer to the surface than Zc for 
Al(lll) (as found for the jellium zq and Zc), and is closer 
than Zq for jellium (as for the Al(lll) and jellium Zc). 
The first trend can be thought of as resulting from the 
difference between the XC hole and the screening charge 
density caused by a point charge, and has been discussed 
for jellium by Eguiluz and Hanke in terms of 'elec- 
tron overlap effects'. The shift inwards for Al is how- 
ever somewhat less than for jellium. The second trend 
is caused by the effects of atomic structure. Whereas 
the jellium model predicts that zq and Zc should be a 
fixed distance from the geometric edge, Zc is seen to be 
more closely tied to the position of the outermost layer, 
with the screening charge density centred just above the 
atoms. The value of zq for Al(lll) suggests that the XC 
hole is similarly modified from the cylindrically symmet- 
ric case of jellium p2[ , (although the size of this shift is 
also somewhat smaller than that for Zc) and that zq will 
thus also tend to follow the position of the outermost 



crystal plane. 

The issue of the physical origin of Vim has provoked 
considerable controversy. Eguiluz et al. |^ found that 
the local potential at a jellium surface contained an ex- 
change part with a dependence, and thus also as- 
cribed the image tail to correlation, in contrast to earlier 
work by Harbola and Sahni [ p3[ , who obtained the image 
potential in Vxc as the work done against a bare exchange 
hole. This discrepancy was explained more recently by 
Solomatin and Sahni who showed analytically that 
the asymptotic exchange image form is given by a highly 
delocalised exchange hole in the semi-infinite metal, and 
predicted the observed faster decay of the exchange part 
in the case of a slab geometry. Thus, although technically 
exchange alone can provide the image tail for a genuinely 
semi-infinite metal, inclusion of exchange and correlation 
effects are clearly necessary to provide a correct physical 
(slab-width independent) description for electrons in the 
near-surface region. 

In any case, Vxc is felt only by electrons in occupied 
states, whereas our potential is that felt by excited states. 
The exchange part of our effective local potential shows 
exponential decay (as suggested by Almbladh and von 
Barth [^), rather than power law behaviour. This is 
because the exchange part of the self-energy is given by 
Sx = iGV where G is evaluated at an infinitesimal pos- 
itive time, and thus contains only occupied states. As a 
result (neglecting constant factors) 

, , /rfr'EV'oec(r)C,,(r')^*Qp(r') 
Viocir) oc — . (6) 

Moving into the vacuum, the exponential decay of the 
occupied states in G means that VJ^^ must decay expo- 
nentially unless the denominator falls equally quickly. In- 
clusion of correlation is therefore essential to describe the 
image potential felt by unoccupied electronic states. 

In conclusion, we have evaluated the electronic self- 
energy at an Al(lll) surface within the GW approxima- 
tion. Calculation of the full energy-dependent self-energy 
and solution of the quasiparticle equations are made pos- 
sible by use of the real-space imaginary-time method. 
The resulting effective local potential shows that correla- 
tion must be included to obtain the image potential for 
excited-state electrons such as those involved in photoe- 
mission or low-energy electron diffraction. The image- 
plane position for the many-electron system is closer to 
the surface than that for classical response to external 
fields or charges, and is also significantly modified by the 
atomic structure of the surface. 
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TABLE I. Image-plane positions for Al(lll) surface in 
a.u., relative to the geometric edge, and for jellium with 
rs = 2.07. Results for dynamic electrons (jellium calculations 
of Ref. [9], Al(lll) this work) arc compared with calculations 
based on response of the ground-state density by Lang and 
Kohn [6], Lam and Needs [19], Finnis [8], and Serena, Soler 
and Garcia [24,20] (see text). 
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FIG. 1. Surface state QP wavefunction (full line), obtained 
by the iterative solution of the energy dependent QP equation 
at r, has weight transferred into the vacuum relative to the 
corresponding KS-LDA eigenfunction (dashed line). 
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The XC potential calculated 



from E crosses over to the classical image form in the vacuum 
(best-fit Vim shown). 
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FIG. 3. Exchange and correlation contributions to the 
effective local potential for excited-state electrons at the 
Al(lll) surface. The bare exchange part decays exponentially 
into the vacuum. 
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